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0^ ' The infinite-C/ three-band Hubbard model is considered in order to describe the Cu02 planes 

of the high temperature superconducting cuprates. The charge instabilities are investigated when 
the model is extended with a nearest-neighbor repulsion between holes on copper d and oxygen p 
orbitals and in the presence of a long-range Coulombic repulsion. It is found that a first-order valence 
instability line ending with a critical point is present like in the previously investigated model without 
long-range forces. However, the dominant critical instability is the formation of incommensurate 
charge-density-waves, which always occur before the valence-instability critical point is reached. An 
effective singular attraction arises in the proximity of the charge-density wave instability, accounting 
for both a strong pairing mechanism and for the anomalous normal state properties. 
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^ ■ I. INTRODUCTION 

, The crucial structural elements, present in all the superconducting cuprates, are the Cu02 planes. At half-filling, 

■ i.e. one hole per unit cell in the Cu02 planes, these materials are antiferromagnetic charge-transfer insulators and 
I , the holes mostly occupy the copper sites. The holes added by doping reside on the Cu02 planes and have a large 

' ^ amplitude on oxigen, thus showing that the strong hole-hole repulsion on copper is a relevant feature. Upon doping 
g , the system is driven toward a paramagnetic metallic phase which becomes superconducting at low temperature. The 
^ • metallic phase above Tc presents many anomalous features, which contrast with the usual behavior of the normal 

'—'I Fermi liquids. 

^ • The failure of Landau Fermi-liquid theory in the metallic phase of the cuprates has been ascribed to singular 
^ interactions arising in the proximity of some critical point at zero temperature (quantum critical point, QCP) [0-^. 
^\ ' More specifically the complex behavior of these systems was recently interpreted in terms of a proximity to an 
incommensurate charge-density-wave (ICDW) transition located at zero temperature near the optimal doping. 
In the quantum critical region above it no energy scale besides temperature rules the physics and strong critical 
fluctuations are responsible for both non-Fermi liquid behavior and large superconducting critical temperatures. As 
• soon as superconductivity takes place the CDW instability is hindered and can only be recovered by destroying 
^\ \ the superconducting coherent state, like, e.g., in transport experiments under strong magnetic fields 1^. In the 
^ ' underdoped compounds, instead the instability would occur at finite temperature were it not for the quenching due to 
\ superconducting local fluctuations which can give rise to the appearance of charge and spin gaps of d-wave symmetry 
^ ■ as experimentally found ||lO|-[T^ at a temperature T* above T^. According to this proposal 0J^, the underlying 
\ \ hindered charge instability provides the temperature dependent pairing potential needed to explain the high crossover 
^ ' temperature T* of the gap formation and the peculiar doping dependence of T* (strongly increasing with decreasing 
Q \ doping, JlO|-|T^) with respect to the value of the charge gap at T=0 (nearly constant with decreasing doping, |ll|,|l^). 
O ' The existence of a QCP near optimal doping is supported by several experimental findings. In particular recent 
transport measurements in La2-xSrxCu04 (LSCO) under strong magnetic field ||^ investigated the normal phase 
of these systems when superconductivity is suppressed. This analysis shows the existance of a QCP near optimal 
. , doping. Indications in this sense are also provided by neutron scattering [l^ and by the qualitative changes of 

■ behavior at optimal doping detected by optical spectroscopy NMR [|^, susceptibility neutron scattering 
I p^ , photoemission ]lO|-p^, specific heat thermoelectric power jl^. Hall coefficient resistivity |9 p^ , |2l[ |. It 
is also suggestive that several quantities (resistivity. Hall number, uniform susceptibility) display a scaling behavior 
with a typical energy scale, which vanishes at optimal doping p^-p^. 

Many indications exist that the above QCP involves charge ordering. Generically, the sizable doping at which the 
QCP occurs suggests that charge degrees of freedom are substantially involved in the ordering phenomenon. A direct 
observation of charge-driven ordering was possible by neutron scattering |^5|-^7|], in Lai.48Ndo.4Sro.i2Cu04 where the 
related Bragg peaks were detected. For this specific compound the low temperature tetragonal lattice structure pins 
the CDW and gives static order and semiconducting behavior (see also the case of Lai.88Bao.i2Cu04). Increasing the 
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Sr content at fixed Nd concentration, the pinning effect is reduced leading to metallic and superconducting behavior. 
In this latter case, the existence of dynamical ICDW fluctuations is suggested by the presence of dynamical incommen- 
surate spin scattering, although the charge peaks are too weak to be observed. In this regard, also the La2-xSrxCu04 
is expected to display dynamical charge fluctuations with doping-dependent spatial modulation as indeed observed 
in the magnetic scattering |28|. ICDW have been proposed from extended X-ray absorption fine structure (EXAFS) 
experiments both in optimally doped LSCO a nd Bi2Sr2CaCu208+x (Bi-2212) Q. Superstructures have also 
been detected in Bi-2212 from X-ray diffraction 

The occurrence of an ICDW-QCP is theoretically substantiated by explicit findings within the single-band infinite- 
U Hubbard model in the presence of an Holstein electron-phonon coupling. In the absence of long-range (LR) 
coulombic forces and for sizable but realistic electron-phonon coupling, this model displays phase separation (PS) 
into macroscopically large half-filled insulating regions and metallic hole-rich doped regions. When the LR Coulomb 
repulsion between holes is included the q — charge instability is prevented. Nevertheless a phase-separation tendency 
remains in the system and shows up on a local basis: The PS region of the phase diagram is replaced at small doping 
by a smaller ICDW instability region characterized by a non- vanishing incommensurate wavevector qc- 

As discussed above, the occurrence of an incommensurate charge instability and of the related singular attractive 
interactions is a key point, which could provide a unified explanation for both the normal and superconducting 
properties of the cuprates. In this regard it seems quite important to investigate how generic and robust is the above 
scenario based on an ICDW instability arising from the competition between PS and LR repulsion. 

As far as PS is concerned, this is a common feature of strongly correlated electron systems: It has been found in 
models with nearest- neighbour Coulomb interactions |32-3^, in magnetic models [^ , ^ , and in models with phononic 
interactions . Indeed a strong on-site correlation drastically renormalizes the kinetic energy, which would tend 

to delocalize the carriers into Bloch quasiparticle states. Then short-range interactions (magnetic, phononic, nearest- 
neighbor coulombic.) introducing effective attractions between the carriers may dominate and give rise to charge 
aggregation in highly doped metallic regions together with charge depletion in spatially separated (ordered) regions 
with no itinerant charges. 

As pointed out by Emery and Kivelson |^ LR Coulomb forces effectively oppose the separation of charged particles 
suppressing long-wavelenght density fluctuations. This may lead to either dynamical slow density fluctuations ^] 
or static ICDW Q, Q, [|^, However, whereas the occurrence of PS is well established in many models of 

strongly interacting electrons, the occurrence of ICDW from LR coulombic repulsion frustrating PS in a microscopic 
electron model was only discussed in a single-band Hubbard-Holstein model and only marginally considered 

in a three-band extended Hubbard model |Q. Then, although the physical ingredients seem fully general, to gain 
further understanding of the ICDW-QCP scenario, it is desirable to investigate other theoretical models of strongly 
interacting electrons in the presence of LR forces. In this way one can both (a) explicitly check the generic character 
of the physical mechanism and (b) highlight the specific features arising from the various interactions. In this regard, 
the simple single-band Hubbard-Holstein model having already been considered, the natural extension seems to be a 
multiband Hubbard model with purely electronic interactions. 

The modelization of Cu02 planes by the three-band Hubbard model, in which the copper and oxigen degrees 
of freedom are explicitly taken into account, was proposed by Varma et al. together with the inclusion of the 
nearest- neighbor Coulomb repulsion in order to stress the role of the charge-transfer (copper-oxigen) fluctuations . 
Littlewood et al. , , using a weak-coupling approach (for a finite not too large repulsion on copper site) have 
found that the presence of a nearest-neighbour Coulomb repulsion V of order of the bandwidth can drive the system 
toward a valence instability (VI), characterized by strong charge fluctuations between the copper d orbitals and the 
oxigen p orbitals. If the total number of particles is allowed to fluctuate, the VI has been shown t o be embedded 
in an unstable region where the system phase-separates. Also in a strong-coupling approach ]3^ , ||34| the mean- field 
solution at fixed number of particle displays a VI line ending at the so-called VI point, where the charge-transfer 
excitonic mode completely softens. However, when the number of particles can fluctuate, this line is surrounded by a 
region of negative compressibility, i.e. by a region in which the system phase-separates. Therefore PS always occurs 
both in weak and strong coupling before any VI takes place. 

In this paper we address the problem of introducing the LR forces in a three-band Hubbard model. In this way we 
study how the charge-transfer fluctuations, which are able to provide a strong attraction in the particle-hole channel 
and lead to PS in the pure short-range case, are affected by a Coulomb potential. In section II we introduce the 
model. In particular we work within a strong coupling approach: An infinite repulsion between holes on the same 
copper site is handled by means of a standard slave-boson technique. Fluctuations around mean-field values are taken 
into account up to the gaussian order using a systematic expansion. This section is rather technical and the 
uninterested reader can skip it only retaining the form of the model provided by the pure short range terms [Eq. (Q) 
with Ud = oo] and by the LR potential [Eqs. (^-(||)]. 
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In section III we will show that, as in the LR Hubbard-Holstein model, the phase diagram of three-band Hubbard 
model in the presence of LR interactions displays a finite-g instability related to CDW formation. We find a singular 
scattering between quasiparticles near the line on which the generalized density-density correlation functions diverge 
at a finite qc- By contrast in the proximity of the VI line no singular scattering is found. Only at the isolated 
VI point the interaction takes a singular form due to vanishing of the excitonic energy, but this has no physically 
relevant consequences since the VI point is always inside the unstable region with respect to the ICDW formation. 
This conclusion differs from the results obtained by Varma in a three-band Hubbard model , where an instability 
related to intracell charge-transfer currents was found, together with singular interactions. However, in the model 
considered by Varma a relevant role was played by general symmetry properties, which are absent in our case, thus 
rendering a direct comparison between the two models and results questionable. In section IV wc discuss the effective 
interaction between quasiparticles, whereas in section V we draw our conclusions. 



II. THE MODEL 
A. Introduction of LR potential 

The Hamiltonian of the three-band Hubbard model in the presence of a local repulsion on the copper site Ud and 
a nearest-neighbour repulsion V reads 

ia iaoL—x.y iaiq—^x^y 

-ippYlp\+xAP^-y<r-V^+ya + {i^i + l) + h.c]+Ud^ni^nf^+V nt"!*^' (1) 

ia i icrcr' rj—^x ,y 

Where {die) creates (destroys) a hole in the copper site i, p\^^ (piaa) in the oxigen site i + a, nf^ — dj^dicr, 
nl^ = plrf^Pirja and and are the bare energy levels. We work in hole representation and we denote by 6 the 
number of holes per unit cell added by doping: ^ [i^fa) + i^ia)] = (1 + '^)- 

In the case oi Ud — oo (strong-coupling) no double occupancy arises on copper sites d\^dia < 1. We handle this 
constraint within the standard slave-boson technique |46| - [50| by introducing a new boson degree of freedom labeling 
the empty site df^ dl^bi, di^ b\dia- Since the site can only be either singly occupied by a fermion or singly 
occupied by a boson, the constraint becomes d\^dia + h\hi = 1. In order to avoid perturbative assumptions on the 
coupling constants we use the large- iV expansion, assuming the spin index to run from 1 to A'^ and letting N oo. 
So the constraint is relaxed into 

E<U. + &l6^ = y, (2) 

(7 

showing that bi ^ ^/N. Consistently we rescale the coupling tpd tpd/VN, V — > V/N. Then the partition function 
of the system can be written as a functional integral 



Z = J VdlVd^Vpl^Vpo.a'Db'^'DbVXDXVYeyi^ drS^ , 



(3) 



with 



^ = E 4^ + E ^l.%r + E ^li^ + ^ E A. {^^^ - <^^^) + ^ E i^^ + + n, (4) 



laa—x.y 



1(77]— :Lx,y 



^ E ^dniv) [hd\„Pir^a + h.c)j - tpp^Pl+xa [Pi-ya - Pt+ya + (i ^ i + 1) + h.c] . (5) 
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A is a Lagrange multiplier introduced to enforce the constraint in Eq. (Q). Writing the nearest-neighbour Coulomb 
interaction as 



X)icrr(=±2; y '^Zr ' ~ Sicr "'fcr ' have decoupled this interaction by means of an Hubbard-Stratanovich 
transformation introducing the fields and 1^. Note that in the leading order of the large- iV expansion only the 
Hartree factorization is present, the Fock term being suppressed by a factor due to the sum over the spin. In order 
to consider the Fock factorization also, one needs to add the term 

-^^2 = ^ 4adta'pl,j„,Ptr,a- (6) 

iaa' ri—±x .y 

This term can be decoupled by a complex Hubbard-Stratanovich field Z^,,. However we shall see in Appendix B that 
the Fock term does not modify the properties of the model from a qualitative point of view. So we discard it from 
now on. 

Now we extend the formalism in order to consider the Coulomb LR potential by taking into account the symmetries 
of the underling lattice, that is of the square lattice of copper atoms with a two-oxygen basis along the directions 
(a/2, 0) and (0, a/2). Then we add to the hamiltonian the term: 

^-'^^E"?^*"''(9X, (7) 

qaa' 

where a,b = d,x,y and the sum over a and h is understood, 
reads as 

1 cos (^) cos (^' 

^-\q)^VM\ cos{^) 1 cos(^)cos'(^) ) , (8) 

and Vc is given by 



cos(^) cos(^)cos(if^) 1 



Vc{q) - ^ „ , (9) 



\j {tT^ [cos(afe) + cos{aqy) - 2] - l| - 1 



We define e|| and ex as the dielectric constant in and out the Cu-0 plane, a and 6 the in-plane and out-plane lattice 
constants; V measures the strength of the Coulomb potential. 

Vc{q) is obtained through the solution of the 3-D Laplace equation for a point-like charge in a cubic lattice, projected 
onto the z = plane |39j. The matrix form of $'^''((7) reproduces the structure of the Fourier-space Laplacian operator 
in the orbital space. In the small q limit Vc{q) has the right behavior 

lim Vc{q) ^ -. 

q 

corresponding to a two dimensional plane embedded in a three dimensional space. 
The partition function including the long-range term now reads 



j VdlVd^Vpl^Vpaa'Db''VbVXDXVYVn''exp (- drSj , (10) 



5 = 5+-5]l]^<i>,-i(q)17^ + z^l]X- (11) 

q qa 

where is a real vector field, introduced to decouple the LR Coulombic interaction. Notice that in Eq. (^ we use 
the long-wavelength uniform values for the dielectric constants in order to obtain the repulsive Coulomb potential at 
long-distance. We have also included in the Hamiltonian the explicit n-n Coulomb repulsion V . This term is needed 
to describe the reduced screening present at short distances between localized charges, which leads to an enhanced 
repulsion between n-n sites. 
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B. Large- A'^ Expansion 



In this section we present the formahsm needed to handle the gaussian fluctuations of the boson fields around the 
saddle-point solutions. From now on we choose the radial gauge Q. In this gauge the phase of bi is gauged away and 
only the modulus is considered, while Ai acquires a time dependence. The Hamiltonian of the coupled fermions and 
bosons is then written in a compact form, by introducing a seven-component field A'^ = (i5r, JA, JY, $7") defined 
as the fluctuating part of the fields around the saddle-point solutions 

n = ro(l + Sri) 
\i = -iXq + SXi 
Xi = Xq + SXi 
Y, = -lYo + SY, 

We have not included the (infinite) zero-momentum component of the Coulombic field $1° since we are assuming that 
it is cancelled by the contribution of a uniform charged ionic background. 
The Hamiltonian reads 

H ~ Hmf + Ti-bos + Ti-int, (12) 
where, using the notation '^kaa — idka-,ipxk<nipyka-), T^mf is the mean-field fermion hamiltonian (in units of a = 1) 

/ ed ~2rotpdSm(!f) -2roip<i sin (^^^ ^ 

-2rotpdSm(!^) €p -2tppPk (13) 

^-2rotpdSm(^] -2tppf3k £p j 



Ti-MFik) = 



where Cp = Cp — Xq + Yq and + Aq + + Yq are the renormalized energy oxigen and copper levels respectively, 

and f3k = 2 sin {kr^ / 2) sin {ky/ 2). 

The quasiparticle basis ^'fc is obtained from the unitary transformation, 'ffc — U{k)'^k which diagonalizes the 
mean- field Hamiltonian Hmf, so that 

The boson-fermion interaction term can be written in the form 

kqa kqcr 



The 3x3 boson-fermion interaction vertices A^^{k, q) in the orbital operator basis can be obtained from Eq. ( |ll| ) and 
are shown in Appendix A. The quasiparticle vertices A''''(fc, q) are given by 

A^{k, q) = u[k+^) A^(fc, g)[/t (fc - I) . (15) 

Finally the purely bosonic part of the hamiltonian is 

Hbos ^NY,A>^{q)B>^''{q)A''{-q), (16) 

where the 7x7 B{q) matrix is given by: 

and the only nonzero elements of Bgr are BH = toAq, BI^ = BH = ir§ and B^f = B^,^ — \/2V. 
At leading order in 1/A^, the dressed propagator of the A^ field is 
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D^^^{q,u;) = {A''{q,u^)A%~q, = ^ [2B{q) + n(g,c.)];,^ , (18) 

where Il.^^{q^u)) are the bare polarization bubbles: 

W!^/ A fa{k + q/2) - ffi{k - q/2) ~ ^iv n \ ^^n\ 

^ ('^'"^ - 2-^^ E^{k + q/2) E,[k - q/2) - 

The factor 2 multiplying the bare boson propagator is due to the fact that in the radial gauge the bosonic fields are 
real. Correspondingly, at leading order in the density-density correlation functions are given by 

where 

X%{q,^) = -^^{naa{q)ni3a'{-q))o (21) 

aa' 

are the bare density-density correlation functions, and 

P^piq,^) = ^E<'^"-(9)E*L'^A!;,(fc,g)vI/,.+,^,,)o, (22) 

aa' k^S 

with a = d,Px,Py and /i running on the boson indices. Linearly combining Xa/?) we calculate the total density 
correlation function 

Xnn{q,t^) ^ {{np + iid){np+nd)) (23) 

and the charge-transfer correlation function 

XCT{q,i^) = {{rip - nd)inp - Ud)). (24) 

Within this formalism it is possible to calculate the residual (order 1/A^) effective scattering amplitude between 
quasiparticles in the lowest band (where the Fermi level lays) 

r{k,k',u;) = -A'^,{k',-q)D^^''{q,cjWn{k,q). (25) 

We note in passing that, being the bare bubble defined in Eq. ( p^ ) and the vertices defined in Eq. ( p^ non-singular 
functions of q, in the present RPA-like approximation the effective interaction Eq. ( |2^ ) and the correlation functions 
Eq. ( pO[ ) always diverge together, since any singularity may only arise from the boson propagator entering both 
quantities 



III. PHASE DIAGRAM AND CORRELATION FUNCTIONS 



The three-band extended Hubbard model with short-range interactions has been widely studied in both weak- 
coupling ||4^- ^ , ^ and in semianalytic js^jsj,^ or fully numeric [pT|-^ strong-coupling approaches. Before pre- 
senting the effects of LR forces on this model, we recall some of the most important results for the pure short-range 
case. In particular we summarize the phase-diagram in the Aq — 5 space, Aq = Cp — being the bare difference 
between the copper and oxygen atomic levels, for various given values of the nearest-neighbor Cu-0 repulsion V in the 
more realistic strong-coupling approach. At half-filling a metal-charge-transfer-insulator (MCTI) transition is present, 
for Aq larger than a critical value (Aq + V> 3.3tpd when tpp — 0). Away from half- filling, increasing V the system 
shows phase-separation above a critical value V* {V* ~ 1.63tpd for tpp = 0). If is increased further, two solutions 
with different occupations of the copper and oxigen orbitals can be found solving the self-consistent equations at 
fixed number of particles. Correspondingly to these solutions, the phase diagram, at finite doping, would display a 
"first-order" line between a d-like metal and a p-like metal ending in a "second-order" valence instability (VI) point 
characterized by the softening of the energy of the p — d exciton mode at g = 0. However, if the number of particles is 
allowed to fluctuate the VI line is always inside a region of negative compressibility where the system phase-separates. 
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Then the VI point is not physically attainable. The region of negative compressibility is delimited by a spinodal line 
on which the compressibility diverges. Together with the compressibility, all the static charge correlation functions 
diverge because of the mixing between the charge fluctuations. Along the spinodal line, the exciton mode remains 
massive (i.e. uJexc{q = 0) ^ 0), even if the charge-transfer (CT) correlation function is divergent. The instability is 
instead accompanied by the overdamped zero sound acquiring a vanishing and then negative velocity. These results 
are only quantitatively, and not qualitatively, changed if the Fock term in Eq. is considered . 

When the LR potential of Eq. (|^) is added to the short-range extended Hubbard model of Eq. (||) , it provides a 
huge electrostatic cost to the long-wavelength charge fluctuations thereby preventing PS. Nevertheless, an instability 
region characterized by the divergence of the density-density correlation functions at a finite Qc is still possible and is 
indeed present, as in the single-band Hubbard-Holstein model, where an ICDW stripe instability was found |,||. In 
Fig. 1 and (2) we show two phase diagrams for two different values of the LR Coulomb force strength V: V = l.Gtpd 
and V = 8tpd with tpp = 0.2tpd and V — 2.3tpd These values of V would correspond to a repulsion between holes 
on two nearest-neighbor copper sites of the order O.ltpd and 0.5^^^) respectively. The phase diagrams are determined 
by identifying the divergences of the various correlation functions within the expansion at the leading order . 
The finite-g instability is indicated by the dashed line. When the strength of the Coulomb potential V is larger, as in 
Fig. 2, the unstable region shrinks around the VI point. The dynamical properties of the extended Hubbard model 
at g = are not affected by the introduction of the LR forces. A remarkable consequence is that the VI point, which 
can be characterized by xcriq = 0,(jJ ^ 0) = oo occurs for the same Aq and S of the purely short-range case. In 
the presence of a (large) V, part of the first-order VI line lays outside the region of the finite-g instability. However 
we have checked that the system does not show any critical behavior along this line due to the fact the energy of 
the exciton-mode remains finite on it and no divergences arise in the response functions as well as in the effective 
scattering amplitude between quasiparticles. 

The static (w = 0) charge-transfer correlation function xct(9,w — 0) is reported in Fig. 3 outside and inside the 
CDW-unstable region, for V = l.Gtpd- Near the instability, a huge peak at g ~ develops in the charge-transfer as 
well as in all other correlation functions and they eventually diverge when the charge-instability line is reached. 

Although it may be influenced by some specific features of the band structure (proximity to a Van Hove singularity, 
partial nesting of some portion of the Fermi surface and so on), the critical qc is not determined by any nesting 
property of the Fermi surface. It rather results from the balance between the tendency towards PS of the short-range 
model and the electrostatic cost imposed by LR forces to charge segregation. Then its modulus depends both on 
doping and, less sensitively, on the Coulomb potential through V. It increases ranging from 0.5 to 0.8 (in units of the 
inverse of lattice constant a) with doping and, slightly, with V. This behavior reproduces the one already found in 
the single-band Hubbard-Holstein model, although in this latter case the value of the modulus of qc was about 20% 
larger. However, this feature is obviously non-universal and is a specific outcome of the band parameters considered 
here. We found that the correlation functions are rather isotropic in the {qx,qy) space. However a detailed analysis 
of the direction along which the instability first occurs shows that the (1, 0) and (0, 1) directions are more favorable 
than the diagonal (1,1) one. 

To clarify the physics underlying the ICDW instability, we studied the behavior of the collective modes at finite q 
and w in the stable region. To this purpose, we analyzed the density-density correlation functions [Eqs. (|3|) and (^]. 
The frequency and momentum dependence of the poles (or resonanaces) of the correlation functions as a function of 
q give indeed the dispersion of the collective modes of the model. Two modes can be identified: The plasmon mode 
and the charge-transfer exciton. Having introduced a LR term proportional to \/q (for small q) the energy of the 
plasmon mode vanishes at g = as y/g, whereas the frequency of the exciton mode at g = stays finite all over 
the stable region of the phase diagram. Due to the mixing, both poles are present in all the correlation functions, 
although with different spectral weight. Whereas the peak of the plasmon is more pronounced in the total density 
correlation function, the excitonic mode has a larger weight in the CT correlation function. 

In Fig. 4 we compare the dispersion of the 2_D-plasmon near and far away the CDW instability. Fig. 4 clearly shows 
that the 2£)-plasmon mode softens at the same g for which the instability takes place. On the contrary, the exciton 
mode remains at higher energies and the CDW instability does not affect much its dispersion. We have also analyzed 
the plasmon mode as a function of the direction of the momentum g in the Brillouin zone. We find that, up to g ~ 0.5 
(in units of the inverse lattice constant), the dispersion is nearly isotropic and depends very weakly on the direction 
of g. In particular, when the system is far from the instability, the dispersion on the (1,0) and (1,1) directions is linear 
(apart from the small region around g = where the q^^^ behavior is found) with a 10% higher slope in the (1,1) 
direction. This analysis illustrates that the occurrence of an ICDW instability in the three-band extended Hubbard 
model is accompanied by the softening of the plasmon mode due to the attractive effective interactions mediated by 
the charge-transfer fluctuations in the presence of a sizable nearest-neighbor repulsion V . The microscopic dynamics of 
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the instability is therefore substantially different from the analogous instability taking place in the Hubbard-Holstein 
model |3^]. In this latter case the electronic plasma mode mixed with an optic phonon mode. The finite momentum 
instability was then accompanied by the softening of the phonon mode, which was also responsible for the effective 
attraction eventually driving the instability. However, despite the physically different microscopic driving forces and 
the different dynamical evolution of the collective modes, the similar occurrence of an ICDW instability driven by the 
interplay of PS tendency and LR forces clearly illustrates the full generality of this mechanism. 



IV. EFFECTIVE INTERACTION 



Our analysis now proceeds with the investigation of the effective interaction between quasiparticles on the Fermi 
surface as defined by 

rikp,k'p,q; K{k'F,-q)D^'' {q, u)K'i^{kF,q). (26) 

Fig. 5 displays T{q,uj = 0) for V = OMpd and 5 — 0.2, for two different values of Aq close to the CDW instability, 
as a function of [q^ , 0) . As in the case of the short-range model p3|-|35t , near the charge instability, a large attractive 
interaction between quasiparticles is generated at = 0. However, while in the short-range model the interaction 
diverges negatively at g = 0, in the presence of LR forces, T{q,uj = 0) stays finite at g = and diverge negatively at 
(7c, m correspondence to the divergence of the correlation functions. This is rather natural since at leading order in 
the singular behavior of both quantities is related to the singular behavior of the boson propagator. 
In the proximity of the instability, T{q,Lo) can be fitted by the following expression 

r(g, Lo) = U- : . (27) 

^ ' {q - qcf + -ljuj ^ ' 

U describes the residual, nearly fc- independent repulsion mediated by the slave bosons: within the large-N slave-boson 
formalism, the infinite U repulsion between the bare fermions, is reduced to a rather weak residual repulsion between 
the Fermi-liquid quasiparticles. The imaginary term in the denominator of Eq. (|2^) is proportional to cu through 
the damping coefficient 7 and reproduces the behavior of the imaginary part of the bare density-density polarization 
bubble of the quasiparticles /TOX°„(w,q) at small frequencies and finite momenta close to qc- This means that. 



despite the complicated structure of the scattering amplitude Eq. (26), in the proximity of the charge instability, the 
interaction between quasiparticles has a simple RPA-like form 

r(g, cj) = u- . ^ ^ , , — r (28) 
1 + r,^(g,w)x^„(g,w) 

with Tt^(q,uj) being an effective dynamical interaction between the quasiparticles. Tt^(q,uj) is dominantly attractive 
at q ~ and is mediated by virtual high-energy processes (mostly interband transitions). The behavior of as 
a function of 6 — 6c is plotted in Fig. 6. It vanishes as a{S — 6c)'^'^ with 2z/ = 1. On general grounds, this is 
what one expects approaching a Gaussian QCP. Since the leading-order expansion bears resemblance with a 
RPA resummation, it is quite natural that the propagator of the critical fluctuations assumes the universal form 
of Eq. (|2^), which is also found in the context of the Hubbard-Holstein model |^,^ and in the proximity of an 
antiferromagnetic QCP where the instability also occurs for a finite value of Qaf = {t^, tt). Of course, the specific 
form of Eq. (|2^) could depend from our approximate (nearly mean-field) treatment. Nevertheless the singular nature 
of the interactions mediated by critical fiuctuations is a sound generic consequence. Many physical consequences 
stem from the presence of a QCP related to an ICDW instability, which have already been generically explored in 
previous works [^|j6j. In particular, the singular interactions were related to an anomalously large decay ratio for the 
quasiparticles at the "hot" spots [1/t ~ y^max[ek, T]] , to linear-in-T resistivity in 2D and p ~ T^^'^ in 3D ||] and a 
specific doping dependence of the superconducting critical temperature Tc was found . Similar results are obviously 
obtained, since the above analyses are quite general in so far they are based on the general theory of QCP's. On 
the other hand, the peculiar character of the presently considered model obviously enter in the values of the non- 
universal constants. In particular, whereas the temperature dependence of the mass parameter in Eq. (p7|), can 
only be calculated within a finite-temperature analysis, other parameters are accessible within our T = slave-boson 
calculation. Specifically, as seen in Fig. 6, for typical Hamiltonian parameter, the zero-temperature doping dependence 
of the mass = a{S — 6c) is determined by the value of the coefficient a 2 and turns out to be quite smaller than in 
the phononic Hubbard-Holstein model. As a consequence in the present case the mass term grows more slowly upon 
increasing the doping away from 6c, thus leading to a more extended region with substantial critical fiuctuations. 
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V. DISCUSSION AND CONCLUSIONS 



The scenario here presented in the context of the three-band Hubbard model involves the relevant interplay between 
the plasma and the charge-transfer modes. An obvious completion of our work would then regard the possible 
experimental confirmation of the physical outcomes of the model. Unfortunately, such a connection with experiments 
is not presently attainable on a quantitative level, due to the lack of detailed experimental analyses on collective 
electronic excitations at finite momenta. 

As far as the limitations of our model are concerned, these are manyfold. First of all, we only considered a purely 
two-dimensional electron system embedded in a 3-D inert screening medium, thus neglecting the layered nature of the 
superconducting cuprates. As a consequence, the investigated collective modes only have a purely two-dimensional 
character. It is well known that, as far as the 2D plasma mode is concerned, it qualitatively differs from the collective 
plasma modes of a layered electron gas (LEG) js^] formed by an array of 2D electron gas layers spaced by a distance 
d. In particular, the screening induced by the other layers changes the ^ dispersion of the 2D mode into a continuous 
branch of modes labelled by the wavevector in the direction perpendicular to the layers, with —(ir/d) < Qz < Tr/d. 
Whereas the kz — mode is massive like in 3D isotropic systems (plasma oscillations in phase on all different layers), 
the other modes are acustic with different velocities (w — > oo for kz 0) until at finite momenta in the 2D Brillouin 
zone they all merge in a narrow, infinitely degenerate branch practically indistinguishable from the dispersion of the 
purely 2D mode. However, it is our opinion that some important qualitative indications can already be gained from 
our simplified model, since the LEG modes merge at q ^ ir/d substantially smaller than the typical momenta q of 
the twodimensional Brillouin zone. Therefore, although explicit calculations are still in progress, we believe that the 
physics at finite sizable momenta of the twodimensional Brillouin zone is not very different for the isolated or the 3D 
layered systems and a CDW instability is a natural outcome provided Qc is larger than q. In our case, we found, 
Qc ~ 0.5 (remember that we take the lattice spacing a — 1) thus restricting the validity of our analysis to systems 
where q < 0.5, i.e. d > 6a. 

Other important simplifications of the model here considered are worth being investigated and could prevent a direct 
comparison with real systems. Specifically in the present model, we did not include a small hybridization {t±) between 
adjacent layers, which would render the plasma modes with kz ^ Q (weakly) massive, thus modifing the dynamics at 
small momenta. Furthermore the interplay between the mixed collective modes and the particle-hole excitations was 
considered within a 1/7V expansion bearing a close resemblance with the RPA approximation. Although extensively 
used in the literature, this approximation is not guarantied to provide a quantitative description of the physics of 
the strongly interacting electrons. In particular, at the leading order in both the screening and the damping of 
the collective modes is determined by the quasiparticles, without any account for the incoherent excitations leading 
to incoherent particle-hole continua on much larger energy scales. These excitations could in principle affect the 
dynamics of the system by providing a damping channel for the collective modes even when these latters are outside 
the quasiparticle-quasihole continuum. 

Besides the above limitations of the theoretical model, a comparison with experiments is also made difficult by the 
lack of detailed experimental analyses of electronic excitations at finite momenta. In particular the extensive wealth 
of data provided by neutron scattering experiments could be related to the above discussed physical effects only in the 
case of a strong enough coupling of the CT and plasma collective modes with the lattice or spin degrees of freedom. 
This could render too indirect the access to the physics of the electronic collective modes, with probably too little 
intensity to get detectable effects. 

On the other hand a direct access to electronic excitations at finite momenta is given by electron-energy-loss 
spectroscopy. However, the presently available analyses, do not seem to us conclusive in order to confirm or to 
disprove any anomalous dynamical behavior of the collective plasma modes. In particular, our analysis showed (see 
Section III) that the anomalous softening of the plasma modes at qc occurs in a rather narrow range of doping 
5 ^ 5c ^ Sopt and momenta q ^ qc- This by itself would require a rather detailed and dedicated search. Moreover, 
it is to be remembered that (local) superconducting pairing "quenches" the ICDW instability, thereby preventing a 
complete softening of the plasma mode. Then, presumably the only consequence of the unachieved charge instability 
would be a partial softening of the mode, likely accompanied by an increase of incoherent spectral weight at low 
energies near the critical qc- This latter effect should be separated from the usual substantial background, which 
likely will render its detection quite a difficult task. 

Despite the above theoretical and experimental limitations, the present work provides a definite contribution to the 
theoretical substantiation of the ICDW-QCP scenario. Indeed we proved that LR coulombic forces do not completely 
stabilize the tendency towards charge instability in the three-band Hubbard model extended with nearest-neighbour 
Coulomb repulsion. In particular we explicitely demonstrated the existence of an ICDW instability for a substantial 
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parameter range of the model and for realistic doping values. Moreover, we showed that, although the VI is still 
present in the phase diagram and remarkably the excitonic softening at g = at the VI critical point is not influenced 
by the LR forces, the ICDW instability is dominant as it embeds the VI point thus providing the only possible source 
of critical singular scattering. 

Although we have presented the results in the strong-coupling limit, the qualitative behavior of the model does not 
change in the weak-coupling case (we report the calculations in Appendix B). We have also checked that both in the 
strong- and the weak-coupling approach the inclusion of a Fock term does not change qualitatively the properties of 
the system: The CDW region is still present, surrounding the VI point and no critical scattering is found on the VI 
line. 

These results support the idea that incommensurate CDW is a common feature of correlated-electron systems 
irrespective of the microscopic interaction mechanisms. 

Indeed a similar scenario was also found in the framework of one-band Hubbard model with phononic interactions 
p9| . As far as the collective modes are concerned, their dynamical behavior is obviously related to the specific 
interactions: In the present case the instability is characterized by the softening of the plasma mode, whereas in the 
Hubbard-Holstein model the CDW instability is accompanied by the softening of the phonon mode, leading to an 
instability of the underlying lattice structure. It is quite remarkable that such different dynamics of the collective 
excitations of the systems, reflecting the underlying microscopic differences of the two models, eventually lead to the 
same generic conclusion as far as the existence of the ICDW critical point is concerned and the concomitant presence 
of singular scattering. Indeed, the effective interaction is found in both models to have the same singular behavior 
as a function oi S — 6c and to display an imaginary part linearly dependent on the frequency. In both models, the 
singular effective interaction is found to be attractive in the particle-hole channel for a sizeable range of q around Qc 
in the stable region close to the CDW instability. The scenario presented so far does not include superconductivity. 
However, it is rather obvious (and it was proven directly in the Hubbard-Holstein model [^) that the strong effective 
attraction between quasiparticles, also shows up in the particle-particle channel and leads to (local) Cooper pairing. 
The occurrence of superconductivity (or of superconducting pairs without long-range coherence) would definitely 
modify the above scenario in so far it provides an alternative to ICDW. A true long-range CDW order is actually 
realized overcoming the quenching pairing tendency only when a commensurability condition is realized which pins the 
charge fluctuations. When commensurability effects do not occur, ICDW fluctuations and superconducting pairing 
interplay and compete, with strong effects also on the magnetic response of the system. On the one hand ICDW 
fluctuations create dynamical charge depleted stripes, where magnetic correlations may subsist even near optimal 
doping. On the other hand the (local) Cooper pairing induced by the singular ICDW scattering can be responsible for 
the charge- and spin-gap effects which arise in the underdoped cuprates. This complex interplay of ICDW, Cooper 
pairing and magnetism is not of our concern here, but is surely a most interesting (and difhcult) subject which is 
relevant for the understanding of the cuprates [RJ§] . 
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APPENDIX A: BOSON-FERMION VERTICES 



in Eq. 0: 
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K'^{k,q) = icos(f ) , h.^{k,q) = 0, h.^{k,q) = i , A^(fc,g) = 0. 
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APPENDIX B: TWO-BAND MODEL IN WEAK-COUPLING APPROACH 

In this Appendix, we present an analytic treatment of a simplified model, showing that (i) the energy of the exciton 
mode at g = (denoted by Wexc) is not affected by LR interactions so that the softening of this mode (VI point) 
takes place in the same point of the phase diagram and that (ii) , once the pure SR system is unstable towards phase 
separation, the system with LR forces always displays a finite-g instability for any value of V . In order to keep the 
formal structure more transparent we confine the calculations to the weak-coupling case |60|] . 

The main difference with respect to the strong-coupling case is given by the absence of the slave-bosons Tq and Xq 
needed to treat the infinite on-site repulsion. To further simplify the treatment we also neglect the direct oxigen-oxigen 
overlap {tpp = 0). In this case one combination of oxygen orbitals does not hybridize with the copper d orbitals, and 
gives rise to a flat non-bonding band. Interband processes involving this non-bonding band are decoupled in the 
small-g limit and, even at finite momenta, do not play any qualitatively relevant role. Therefore, for the sake of 
simplicity, we will only consider the simplified two-band Hubbard model. 

The form factor 7^ = \J siv? {kx/2) + sin^ (^j//2) can be introduced to relate the Fourier transform pk of the bonding 
oxygen orbital combination pi, to the Fourier trasform of the copper d orbitals. The resulting mean- field Hamiltonian 
is given by: 



where td and are the atomic levels, including the Hartree shifts. The mean- field hamiltonian in Eq. (Bl) can be 
diagonalized by a standard unitary transformation, giving the two bands E±{k) = \ ep + td± ^A^Tsi^/Tl] /2. We 
define A = Cp — e^. 

We are interested in the analysis of the physical response functions x±,±{<1t'^) for the total density (-I-) and CT (— ) 
fluctuations. A simple analytic progress can only be made in the small-q limit, to which we will confine our treatment 
here. In this limit the complicated internal structure of the LR interaction [see Eqs. (0)-(|9|)] greatly simplifies, and 
the cosine form factors in the nearest-neighbor interaction [see the expressions of the A'^' vertices in Appendix A] 
drop. Moreover, in the spirit of weak-coupling approach, we treat the interaction terms in Eq. ([T]) by means of a 
standard Hartree decoupling. (We will also show below that adding the contribution of the Fock decoupling does not 
change qualitatively the behavior of the model). 

The Hartree decoupling of the n-n coulomb repulsion is the following 

icrcr' ia 

where and are the local values of the copper and oxygen density per spin. 
The resulting interaction matrix (in the total-density and CT basis) reads 

with Vc{q) ~ V/q and where the local U interaction on copper was also introduced. At gaussian level, the response 
functions are given by the matrix relation 

Hn((7,(i;) , 

* ^ = 'y, B4 

l + Uoiq,Lj)Viq) 

involving the bare-bubble matrix: 
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(B5) 



where G°J^p{k,v) is the matrix of the Fermionic Green functions in the orbital basis {a, (3 = d,p) and A"(fc) are the 
vertices, which couple the p and d fermions to the density fluctuations 



A+(fc) 



1 
1 



A-(fc) 



1 

-1 



(B6) 



The value of the CT exciton frequency at zero momentum can be extracted from the pole of the CT correlation 
function, since, due to particle conservation, the total density (+) decouples from the dynamics. This shows up in 
the q Lo ^ Q behavior of the bubbles containing the total-density vertex: 



limn+"(g,c.)~(?2^ 



(B7) 



(a = +, — ). Then, since Vc{q) ^ V/q always couples to a total-density vertex oHIq", the diverging Coulomb potential 
is cancelled by the small-q behavior of the bare bubbles having at least one vertex and does not appear in the 
q = limit of Eq. ( B4 ) . By introducing the notation 



'■'■Ou 



lini n^''(g,a;) 



'■'■Oq 



limn^''(Q~0,w), 



one obtains 



Xct(0,w) = X — (0,^^) 



Bn 



i + nr-y— 



(B8) 



Notice that Vc no longer appears in Eq. (B8), thus showing that ujexciq = 0) 
the pure SR case. 

In the opposite static limit u = and small q 



1/2 

(l + Bf^^y ) is the same as in 



X— (?) = Xcriq) 



n, 



Of; 



{n++no-~(n+-)'}^++ 



V- 



n++n — 



(B9) 
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1 + y—u 



Oq 



Oq 



B 



Oq 



B++n — 

^^Oq ^^Oq 



> 0. 



(BIO) 



Using the relation 



lim B, 



lim 

q^O 



B 



Oq 



(Bii) 



and noting that, for small q's, V^^ is dominated by the LR coulombic term V/q, a simpler form can be obtained for 



Oui 



(B12) 



It worth noting that, due to this relations, the static limit of the correlation function equals the dynamical limit, 
provided ujexc is finite. This fact depends on the presence of the LR forces: in the dynamical limit the intra-band 
processes are switched off and only inter-band transition are allowed; in the static limit both processes are allowed, 
but the LR interactions drastically reduces the intra-band transitions at g = 0, while it allows the inter-band ones. 
This also justifies why the exciton energy does not change introducing LR forces: ujexc is indeed involved in inter-band 
processes which are not affected. 

The zeros of the denominator of Eq. (B12) can be obtained through the equation 
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B'^q^) - A'-q = 0. 



(B13) 



with = Alf/n++ 

instability 



If 



< A^ / BV , Eq. (B13) has two real solutions showing the instaurance of a finite-g 
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exc * 



A^ 
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A^ 
VB^ 



(B14) 



Approaching the VI point upon changing the doping, ujexc decrease. Eventually it reaches the finite critical value uj^^^ 
where the charge-transfer (and any other) correlation function diverges at finite-g and the instability takes place. 
From the above result we see that there is a direct connection between the lowering of the exciton energy Wexc smd 
the developing of a region of instability with the finite-g softening of the plasmon mode. A similar behavior can be 
found in the SR model p3[ |, where the q = instability is driven by the partial softening of the exciton mode which 
pushes the zero-sound into the continuum eventually leading the system to a phase-separation instability. 

Finally, we show that the inclusion of a Fock contribution does not change the qualitative behavior of the model. 
For the sake of simplicity we limit ourself to the weak-co upli ng approach, but the same conclusions hold also in the 
strong-coupling case. In this case we have to add to Eq. (B2) the following terms: 



(B15) 



We introduce two suitable linear combinations of the Fock bosons: 

A^ 
B,= 



Considering a mean-field and fluctuating part for these bosons {Ai = j4o(1 + SAi), Bi = Bo{l + 5Bi)), two new vertex 
are needed in the bare-density-correlation functions (we note in passing that the ^-vertex has the same structure of 
the r-vertex in the strong-coupling approach) : 





-2VAo-/k 



-2VAo^k 




2VBo-fk 
~2VBojk 



(B16) 



Due to the fact that the boson Bi decouples from the others we can easily write the dynamical limit of the correlation 
hmction XCTil) including the Fock terms: 



Xct(w) = X — (^^) 



n 



(B17) 



where Vaa = —2V. This result is the same as the one obtained without the Fock term apart from a contribution 
proportional to Vaa which simply renormalizes the exciton- mode energy cj^^ . 
Furthermore the statical limit is: 



X — (?) = Xcriq) = - 



{n++no--(n+-)'} 



{noV 



V- 



n++n — 



Vaa 



TT + +TTAA 



V++ + Al 



(B18) 



HF 



This expression can be cast in the same form of Eq. ( |B9[ ) provided ujexc is replaced by and so all the conclusions 
we have drawn after Eq. (B12) still hold. 
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[60] The strong-coupling extension of the weak-coupling treatment is rather straightforward, since, besides the total density 
and the CT fluctuations, only one additional bosonic field, r^, is coupled to the fermionic fluctuations, thus giving rise 
to additional response fuctions. All other fields, Ai, Xi, Yi, and Hi, couple to the total density or to the CT density 
fluctuations (or to their combinations) , already considered in the weak-coupling approach presented below. 

Figure Captions 

Fig. 1: Phase diagram in the plane e° — eJJ vs 5 for tpp = Q.2tpd, V = "^-^tpd. The Coulomb interaction strenght is 
taken so that the repulsion between two nearest-neighbor copper atoms, Vc{d— d), is O.ltpd- 

Fig. 2: Phase diagram in the plane e° — vs 5 for tpp = 0.2tpd, V = 2.3tpd. The Coulomb interaction strenght is 
taken so that the repulsion between two nearest-neighbor copper atoms, Vc{d— d), is 0.5^^^. 

Fig. 3: Charge- Transfer correlation function Xcriq, 0) vs momentum for S = 0.12, V = IMpd and e° — = l-835tpd 
(+). - - IMtpd (□), e°p-e° = IMBtpd (x), e° - e° = IMtpd (o). The critical value is - e° ^ l.SMpd. 

Fig. 4: Plasma frequency iVpi vs momentum for S = 0.12 and e° — = l.%2>tpd (o), ^p — ^d^ 1.845^),^ (-I-); the 
critical value is Cp — ~ l-85tpd 

Fig. 5: Effective Interaction T(q,Q) vs momentum for 5 = 0.2 and e° - = imtpd («), Cp - = l-922fpd {+). 
The CDW instability corresponds to — eJi ~ l.'d2Ztpd 

Fig. 6: Doping dependence of the mass parameter at T = near the critical point. 
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